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ABSTRACT 

With the NEOWISE portion of the Wide-field Infrared Survey Explorer 
(WISE) project, we have carried out a highly uniform survey of the near-Earth 
object (NEO) population at thermal infrared wavelengths ranging from 3 to 22 
/im, allowing us to refine estimates of their numbers, sizes, and albedos. The 
NEOWISE survey detected NEOs the same way whether they were previously 
known or not, subject to the availability of ground-based follow-up observations, 
resulting in the discovery of more than 130 new NEOs. The survey's uniform 
sensitivity, observing cadence, and image quality have permitted extrapolation 
of the 428 near-Earth asteroids (NEAs) detected by NEOWISE during the fully 
cryogenic portion of the WISE mission to the larger population. We find that 
there are 981±19 NEAs larger than 1 km and 20,500±3000 NEAs larger than 
100 m. We show that the Spaceguard goal of detecting 90% of all 1 km NEAs 
has been met, and that the cumulative size distribution is best represented by a 
broken power law with a slope of 1.32±0.14 below 1.5 km. This power law slope 
produces ~ 13, 200±1,900 NEAs with D >140 m. Although previous studies 
predict another break in the cumulative size distribution below D ~50-100 m, 
resulting in an increase in the number of NEOs in this size range and smaller, we 
did not detect enough objects to comment on this increase. The overall number 
for the NEA population between 100-1000 m is lower than previous estimates. 
The numbers of near-Earth comets and potentially hazardous NEOs will be the 
subject of future work. 
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Introduction 



The near-Earth objects (NEOs) are a population of interest for a wide range of scientific 
investigations and practical considerations. The NEO population (defined as asteroids or 
comets with perihelion distances q < 1.3 AU; http://neo.jpl.nasa.gov/neo/groups.html) is 
thought to be made up of both asteroids and comets ranging in size from objects tens of 
kilometers in diameter (Shoemaker 1983) down to dust grains. Although nearly ~8,000 



NEOs have been discovered to date at all size ranges {http://neo.jpl.nasa.gov/stats), this 



number represents only a fraction of the total population thought to exist (Bottke et al. 



2002). Insights into several fundamental questions in planetary science can be gained from 



large-scale studies of NEO orbital distributions and physical properties such as diameters 
and albedos. 



1.1. Origins and Evolution of NEOs 



NEOs are a transitory population, with dynamic lifetimes of 10 to 10 years (Morbidelli 



& Gladman 1998). They must therefore be continually replenished in order to maintain the 



population we observe today. It has been shown that NEOs are likely to have been delivered 



(Wetherill 


1988 


Rabinowitz 


1997a|b 


Bottke et al. 


2002) 



such as the z/g secular resonance and 3:1 mean motion resonance with Jupiter. The near- 
Earth asteroids (NEAs) in the NEO population are thought to have been pushed by the 
Yarkovsky thermal drift force from different parts of the Main Belt into resonances that 



may preferentially inject them into the inner Solar System (Bottke et al. 2002). Near-Earth 
comets (NECs) also constitute a part of the NEO population and are thought to originate 
from a variety of regions, including the Kuiper belt, the scattered disk, the Oort cloud, and 
the Jupiter Trojan clouds ( Levison fc Duncan|1997 Duncan fc Levison|1997 Weissman|1996 



Levison et al. 1997). Efforts to more precisely constrain the origins of the NEOs have been 
hampered by the relative paucity of observational data. Surveys such as the Lincoln Near- 
Earth Asteroid Research Program, the Near Earth Asteroid Tracking Program, the Lowell 



Near-Earth Object Survey, the Catalina Sky Survey, PanSTARRS, and Spacewatch ( 


Stokes 


et al. 


2000 


Helin et al. 


1997 


Koehn & Bowell 


2000 


Larson 


2007 


McMillan 


2007 


) have 



identified NEOs, but they observe in visible wavelengths, so they are not able to sense low 
albedo NEOs as effectively, nor can visible light measurements provide strong constraints on 
diameters, albedos, or other physical characteristics. Size and albedo distributions for NEOs 
derived from these visible light surveys remain uncertain (particularly for sub-kilometer 
objects) since these models rely upon the absolute magnitude (H) through an assumed 



albedo to derive size (Stuart & Binzel 2004; Morbidelli et al. 2002; Shoemaker et al. 1990 
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Luu & Jewitt 1998). Studies linking meteorite falls to NEO source regions in the main 



asteroid belt have been performed (e.g. Thomas & Binzel 2010; Vernazza et al. 2008), but 



meteorites have been heavily processed by passage through the Earth's atmosphere, which 
also introduces biases in composition and size. 



Bottke et al. (2002) were limited in their ability to model the NEO population by the 



fact that the Spacewatch data upon which their results were based provided them only a very 
small sample of optically-selected NEOs. The paucity of physical characterization data for 
large numbers of NEOs has complicated efforts to identify weaker source regions. Models of 



NEO origins such as Bottke et al. (2002) can be tested and improved by obtaining a sample 



of NEOs with accurately measured diameters and albedos, well-known orbits, and well- 
understood survey detection biases and efficiency. Infrared observations of all classes of minor 
planets are useful for determining size and albedo distributions, as well as thermophysical 
properties such as thermal inertia, the magnitude of non-gravitational forces, and surface 



roughness (e.g., Bhattacharya et al. 2010 Tedesco et al. 2002 Trilling et al. 2010a Harris 
et al.||2009l). 



1.2. NEOs As Impactors 



Asteroids and comets have impacted Earth throughout its history, profoundly altering 



the course of life on Earth (Alvarez et al. 1980; Hildebrand et al. 1991). Currently available 



estimates of the size distribution of potentially hazardous NEOs lead to the assessment that 
objects capable of causing global catastrophe (larger than 1 km in diameter) are thought 



to impact approximately every 700,000 years (National Research Council Report 2010). 



However, smaller objects can still cause considerable damage (NEO Science Definition Team 



2003 Chapman 2004). Recent simulations suggest that the 1908 Tunguska impact could 



have been caused by an object as small as 30-50 m (National Research Council Report 



2010), which would imply that the potential for damage caused by smaller objects has 



been underestimated. However, this increased risk from smaller impacts is offset by recent 



population estimates (Boslough & Harris 2008; Harris 2008) that indicate that the population 



of asteroids in the size range between several tens to hundreds of meters in diameter may be as 
much as a factor of three less than estimated using a straight-line power law of slope -2.354 as 



proposed by Bottke et al. (2002), or a slope of -1.87. Nevertheless, considerable uncertainty 



about the true NEO size distribution persists due to the lack of accurate diameters, albedos, 
and surveys that are not substantially biased against low albedo objects. Even less is known 
about the subset of NEOs that are considered potentially hazardous. 
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Observations 



WISE is a NASA Medium-class Explorer mission designed to survey the entire sky 
in four infrared wavelengths, 3.4, 4.6, 12 and 22 /im (denoted Wl, W2, W3, and W4 
respectively). In a Sun-synchronous orbit, WISE observed the entire sky in all four filters 
simultaneously using three beamsplitters and a scan mirror. The pre-launch description of 



Liu et al. 


( 


2008 


); I 


Vlainzer 


Wright et al. 


( 


201C 


)). The 



final mission data products are a multi-epoch image atlas and source catalogs that will serve 
as an important legacy for future research. The survey has yielded observations of over 
157,000 minor planets, including Near-Earth Objects (NEOs), Main Belt Asteroids (MBAs), 



comets, Hildas, Trojans, Centaurs, and scattered disk objects (Mainzer et al. 2011a[ hereafter 
MllA). WISE has observed nearly two orders of magnitude more minor planets than its 



predecessor, the Infrared Astronomical Satellite (IRAS; Tedesco et al. 1988 Matson 1986). 
WISE launched on 14 December 2009. The WISE survey began on 14 January, 2010, and 
the mission exhausted its primary tank cryogen on 5 August, 2010. An augmentation to the 
WISE baseline data processing pipeline, "NEOWISE," permitted a search for new moving 
objects to be carried out using the WISE data in near-real time. NEOWISE also allowed 
the individual exposures for each part of the sky (an average of eight on the ecliptic plane, 
but rising into the hundreds near the ecliptic poles) to be archived and made available to 
the public. Exhaustion of the secondary cryogen tank occurred on 1 October, 2010, and the 
survey was continued as the NEOWISE Post-Cryogenic Mission using only bands Wl and 
W2 until 1 February, 2011. 



As described in Wright et al. (2010) and MllA, the NEOWISE survey cadence resulted 
in most minor planets in the WISE sample receiving an average of 10-12 observations over 
~36 hours, although some NEOs were observed dozens or even hundreds of times. NEO- 
WISE detected minor planets using the WISE Moving Object Processing System (WMOPS). 



WMOPS was adapted from the Pan-STARRS Moving Object Processing System (Kubica 



et al. 2007) to account for the different cadence and observing strategy used by WISE. Dur- 
ing the entire period of survey operations, NEOWISE is known to have observed at least 
584 NEOs and is responsible for the discovery of 135 new, previously unknown NEOs. For 
purposes of determining the debiased population characteristics of the NEOs, in this paper 
we consider only those objects detected during the fully cryogenic portion of the survey 
by the WMOPS pipeline using the first-pass WISE data processing pipehne as described 
below. This sample consists of a set of 428 NEOs, of which 314 were recoveries of objects 
discovered by other groups, and 114 were NEOWISE discoveries. By considering this subset 
of NEOs, we can study a population that was observed in a more-or-less uniform fashion. 
This approach facilitates the understanding of the NEOWISE survey biases, allowing us to 
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model the unobserved portion of the NEO population. We will consider the NEOs detected 
during the post-cryogenic portion of the survey in a future work. Figure [T] shows the ob- 
jects detected by NEOWISE using the first-pass data processing pipeline; while most of the 
solar system was surveyed during the fully cryogenic portion of the mission, some of it was 
observed with only bands Wl and W2, leading to a significant drop in sensitivity. For this 
and other reasons, we must debias the survey in order to understand the properties of the 
NEO population and disentangle them from survey biases. 

The observations of the objects listed in Table 1 were retrieved by querying the Minor 
Planet Center's (MFC) observation files to look for all instances of individual WISE detec- 
tions of the desired objects that were reported using the WISE Moving Object Processing 
System (WMOFS). The resulting set of position/time pairs were used as the basis of a query 
of WISE source detections processed with the First Pass version of the WISE data process 
pipeline (Version 3.5) in individual exposures (also known as "Level lb" images) using the 
Infrared Science Archive (IRSA; Cutri et al.||2011 ). In order to ensure that only observations 



of the desired moving object were returned from the query, the search radius was restricted 
to 0.3 arcsec from the position listed in the MFC observation file. A radius of 0.3 arcsec 
was chosen to account for any round-off error between the decimal format used by the WISE 
catalog and the hexigesimal format used by the MFC Additionally, since WISE collected 
a single exposure every 11 seconds, the modified Julian date was required to be within 2 
seconds of the time specified by the MFC. The artifact identification flag cc_fiags was al- 
lowed to be equal to either 0, F or p; this flag indicates that the sources were unlikely to 
have been affected by persistence. We required that the flag ph_qual be equal to A, B, or 
C (this flag indicates that the source is likely to have been a valid detection). As described 
Mainzer et al. (2011b, hereafter MllB) and Cutri et al. ( 2011[ ), we included observations 



m 



with magnitudes close to experimentally-derived saturation limits, but when sources became 
brighter than Wl = 6, W2 = 6, W3 = 4 and W4 = 0, we increased the error bars on these 
points to 0.2 magnitudes and applied a linear correction to W3 (see the WISE Explana- 
tory Supplement for details). Each object had to be observed a minimum of three times 
in at least one WISE band, and to avoid having low- level noise detections and/or cosmic 
rays contaminating our thermal model fits, we required that observations in more than one 
band appear at least 40% of the number of observations found in the band with the largest 
number of observations (usually W2 or W3 for NEOs). The WMOFS system is designed 
to reject inertially fixed objects such as stars and galaxies in bands W3 and WA. However, 
the source density in bands Wl and W2 is ~100 times higher than in bands W3 and W4, 
so it was more likely that asteroid detections were confused with stars or galaxies at these 
wavelengths. In order to remove such confused asteroid detections, we cross-correlated the 
individual Level lb detections with the WISE atlas and daily coadd catalogs. Objects within 
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6.5 arcsec (equivalent to the WISE beam size at bands Wl, W2 and W3) of the asteroid 
position which appeared in the coadded source hsts at least twice and which appeared more 
than 30% of the total number of coverages of a given area of sky were considered to be in- 
ertially fixed sources; these asteroid detections were considered contaminated and were not 
used for thermal fitting. 

During the cryogenic portion of the survey, 231 objects were placed on the MFC's NEO 
Confirmation Page {http://www.minorplanetcenter.org/iau/NEO/ToConfirm.html). Of these, 
115 were designated as NEOs, while the remaining candidates that received visible light 
follow-up were confirmed to be Main Belt asteroids, Trojans, etc. (of the ~300,000 tracklets 
submitted to the MPC, only a handful were rejected as spurious). Follow-up observations 
of NEOWISE-discovered NEOs were carried out by a world-wide network of amateur and 
professional astronomers. Figure [2] lists the numbers of observations of NEOWISE dis- 
coveries contributed by MFC observatory code within 15 days of the first observation by 
NEOWISE. Follow-up observations were critically important for securing good long-term or- 
bits of NEOWISE-discovered NEOs because in general, the arc lengths spanned ~36 hours. 
Of the 115 NEOs discovered during the cryogenic mission, 12 received designations despite 
having no optical follow-up. An additional 22 objects objects were deemed lost because they 
received no optical follow-up, and their arcs were too short to designate. The loss of these 
objects will contribute some uncertainty to the size and albedo distributions computed for 
NEOs from NEOWISE data; this will be discussed further below. 

Figures |3] shows the typical WISE colors for NEOs, Main Belt asteroids and Trojan 
asteroids observed during the cryogenic NEOWISE mission. The NEOs occupy a distinct 
region of color space and can readily be distinguished from other sources such as stars, distant 
galaxies and cool brown dwarfs. Examination of the WISE colors can be used as a means 
to constrain orbits in cases where an object's orbit is uncertain (e.g. because it has only a 



short observational arc). Bhattacharya et al. (2010) suggest that thermal colors can be used 



to identify new candidate minor planets; however, the vast size of the WISE single-exposure 
source lists (which include transient artifacts such as cosmic rays and latent images) makes 
techniques that focus on repeated detections of moving sources the preferred method for 
discovery of new objects. Figure |4] shows W3 — Wi vs. apparent sky-plane velocity and the 
histogram of sky-plane velocities for observed NEOs; most NEOs detected by NEOWISE are 
moving at ~ 0.5°/day. Figure |5] shows the WISE colors compared with heliocentric distance; 
as expected, the colors redden as the objects are located at greater distance from the Sun. 
It may be possible to use the WISE colors and apparent sky-plane velocities to constrain 
orbits in cases where observational arcs are short and orbits are poorly known. A future 
work will explore this further. 



3. Preliminary Thermal Modeling of NEOs 



We have created preliminary thermal models for each NEO observed by WMOPS during 
the fully cryogenic portion of the survey with the First Pass version of the WISE data 
processing pipeline described above; these thermal models will be recomputed when the final 
data processing is completed. As described in MllB, we employ the spherical near-Earth 



asteroid thermal model (NEATM) (Harris 1998). The NEATM model uses the so-called 
beaming parameter t] to account for cases intermediate between zero thermal inertia, the 



Standard Thermal Model (STM) of Lebofsky & Spencer ( 1989 ) and high thermal inertia, the 



Fast Rotating Model (FRM; Lebofsky et al. 1978: Veeder et al. 1989 Lebofsky fc Spencer 



1989). In the STM, r] is set to 0.756 to match the occultation diameters of (1) Ceres and (2) 



Pallas. In the FRM, rj is equal to tt, and the temperature distribution is different than that 
used in NEATM, being independent of longitude. With NEATM, 77 is a free parameter that 
can be fit when two or more infrared bands are available (or with only one infrared band if 
diameter or albedo are known a priori as is the case for objects that have been imaged by 
visiting spacecraft or observed with radar). 

Each object was modeled as a set of triangular facets covering a spherical surface with 



a variable diameter (c.f. Kaasalainen et al. 2004). Although many (if not most) NEOs are 
non-spherical, the WISE observations generally consisted of ~10-12 observations per object 
uniformly distributed over ~36 hours (MllA), so on average, a wide range of rotational 
phases were sampled. Even though this sampling helps to average out the effects of a 
rotating non-spherical object, caution must be exercised when interpreting the meaning of 
an effective diameter in these cases. Figure |6] shows the average peak-to-peak amplitude in 
W3 for NEOs detected by NEOWISE during the fully cryogenic portion of the mission. While 
most objects' amplitudes are ~0.5 mag or less, there is a fraction of the population with 
extremely high amplitudes. These objects may be far from spherical, and so the application 
of a spherical thermal model to these objects may not yield accurate results. These objects 
are discussed in more detail below. 

Thermal models were computed for each WISE measurement, ensuring that the correct 
Sun-observer-object distances were used. The temperature for each facet was computed, 
and the Wright et al. (2010) color corrections were applied to each facet. In addition, we 



adjusted the W3 effective wavelength blueward by 4% from 11.5608/im to 11.0984 /xm, the 
W4 effective wavelength redward by 2.5% from 22.0883 fim to 22.6405 /im, and we included 
the -8% and +4% offsets to the W3 and Wi magnitude zeropoints (respectively) due to 



the red-blue calibrator discrepancy reported by Wright et al. (2010). The emitted thermal 



flux for each facet was calculated using NEATM with the bandcenters and zero points given 



in Wright et al. (2010): the temperature at the anti-subsolar point was set to 3 K, so the 
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facets closest to this point contribute little flux. For NEOs, bands Wl and W2 typically 
contain a mix of reflected sunlight and thermal emission. The flux from reflected sunlight was 
computed for each WISE band as described in MllB using the lAU phase curve correction 
dBowellet al.|1989D . Facets which were illuminated by reflected sunlight and visible to WISE 



were corrected with the Wright et al. (2010) color corrections appropriate for a G2V star. 
In order to compute the fraction of the total luminosity due to reflected sunlight, it was 
necessary to determine the albedo in bands Wl and W2, pjr. This is discussed in greater 
detail below. 

Physical properties were computed for each object by grouping together observations 
with no more than a ten day gap between them. This restriction was imposed to ensure 
that NEOs, which can have significant changes in distance over short times, were modeled 
accurately. 

In general, NEO absolute magnitudes (H) were taken from the MFC's orbital element 
files, and errors on H were taken to be 0.3 magnitudes. However, updated H magnitudes 



were taken from the Lightcurve Database of Warner, Harris & Fravec (2009) for 78 NEOs 



that were observed by WISE. Emissivity, e, was assumed to be 0.9 for all wavelengths (c.f. 



Harris et al. 2009), and G (the slope parameter of the magnitude-phase relationship) was 



set to 0.15±0.1 unless a direct measurement from Warner, Harris & Fravec (2009) was 
available. We note that for some objects with high albedos, the choice of G = 0.15 may not 
be appropriate; this is discussed in more detail below. 

For objects with measurements in at least two or more WISE bands dominated by ther- 
mal emission, the beaming parameter rj was determined using a least squares minimization 
but was constrained to be less than the upper bound set by the FRM case (tt). Figure [T] 
shows the fitted t] values for objects which had measurements in two or more WISE thermal 
bands, along with a best-fitting double Gaussian distribution. The median value of the 313 
NEOs that had fitted 1] was 1.40±0.5. The beaming parameter could not be fitted for objects 
which had observations in only a single WISE thermal band; these objects were assigned 
r] = 1.4 ± 0.5. Figure [T] also shows fitted beaming as a function of phase (a) along with 
the best-fitting linear fit, as well as as a function of heliocentric distance and subsolar 
temperature (T^s). The weighted best-fit relationship between rj and phase a is given by 
T] = (0.00963 ± 0.00015)a + 0.761 ± 0.009. Figure shows that t] is strongly correlated 
with the heliocentric distance at which a given object was observed as shown by the red 
line representing the running median for the NEOs, Mars crossers and Main Belt asteroids; 
the closer an object is to the Sun, the higher its r] value for objects observed near 90° elon- 
gation. Because NEATM assumes that an asteroid's night side contributes essentially no 
thermal flux, objects that are observed at high phase angles (e.g. far away from the subsolar 
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point) may have r] values that vary in model-dependent ways. Since WISE always observed 
near 90° solar elongation, NEOs in particular tended to be observed at high phase angles, 
meaning that part of the night sides were observed. If the nightside temperature is not 
zero as assumed by NEATM, the relationship between t] and heliocentric distance /phase 
angle can vary depending on the real temperature distribution. While the r] parameter can 
serve as an indicator of thermal inertia at low to moderate phase angles, this result suggests 
that attempting to use rj as a gauge of thermal inertia or regolith properties at high phase 
angles must be done with caution, due to the phase-angel dependence of rj. In general, 
detailed thermophysical modeling should be used if possible to unravel the various thermal 
and geometric effects contributing to the observed thermal continuum. 

As described above, bands Wl and W2 consist of a mix of reflected sunlight and thermal 
emission for NEOs, and bands W3 and W4 consist purely of thermal emission. In order to 
properly model the fraction of total emission due to reflected sunlight in each band, it 
was necessary to determine the infrared albedo pm. We have made the assumption that 
Ps.e^m = PiSfim = PiR, although it is known that there are a number of material absorption 



features at these wavelengths (see Gaffey et al. (2002) for an overview; the topic is treated 



in more detail in Mainzer et al. (2011d hereafter MllD)). In order to test the validity of the 



assumption that the infrared albedo is the same in bands Wl and W2, we recomputed all 
the thermal model fits with band Wl only and found that the resulting diameters, visible 
albedos and infrared albedos were identical to within ±10% 1-cr. The geometric albedo pv 
is defined as the ratio of the brightness of an object observed at zero phase angle to that of 
a perfectly diffusing Lambertian disk of the same radius located at the same distance. The 
Bond albedo (A) is related to the visible geometric albedo pv hj A ^ Ay = qpv, where the 
phase integral q is defined such that q = 2 f ^{a)sin{a)da. $ is the phase curve, and q = 1 
for $ = max{0,cos{a)). G is the slope parameter that describes the shape of the phase 



curve in the H — G model of Bowell et al. (1989) that describes the relationship between 
an asteroid's brightness and the solar phase angle. For G = 0.15, q = 0.384. We make 
the assumption that pm obeys the same relationship, although it is possible it varies with 
wavelength, so what we denote here as pm for convenience may not be exactly analogous to 
Pv- We can derive pm directly from the WISE objects that have observed reflected sunlight 
in bands Wl and W2 as well as observations in W3 or W4. Figure |8] shows the ratio of pm 
vs. Pv for the objects that had fitted piR. For the NEOs for which pm could not be fitted, 
we used pm/pv = 1-6 ± 1.0. 



As discussed inMllD and Mainzer et al. (2011e, hereafter MllE), asteroids with dif- 
ferent taxonomic classes can have different pm/pv values despite having similar pv'-, Figure 
[s] shows the comparison between pv and PibJpv for the WMOPS-detected NEOs for which 
Vir/Vv could be fitted. For example, a correlation was found between spectral classes that 
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have neutral or blue visible/near- infrared (VNIR) spectral types, such as C or B types, and 
objects that have lower values of VirIvv-, even though these spectral classes have nearly 
identical py. The figure shows the areas where most of the major taxonomic types group, 
although it should be noted that most of the objects studied in MllD and MllE were Main 
Belt asteroids. The extent to which the same groupings in and pm/pv apply to NEOs 
is yet to be determined. Taxonomic types with red VNIR slopes such as T or D types, 
tend to have higher values of Pm/pv- It is possible that the NEOs shown in Figure [s] with 
low pv and low pm/pv are C or B types, whereas those with high pm/pv and low pv are 
more likely to be T or D types. As discussed in MllD, the current set of spectroscopically 
classified objects is heavily affected by selection biases; adding more objects drawn from the 
NEOWISE sample would improve our understanding of the correspondence between VNIR 
spectroscopic shapes and features with size, pv, and pjr. 



3.1. Thermal Model Error Analysis 

Error bars on D, pv, pm, rj and Tgg were determined for each object by running 50 
Monte Carlo trials that varied the objects' H values by the errors described above and 
the WISE magnitudes by their error bars using Gaussian probability distributions. The 
minimum magnitude error for all WISE measurements was taken to be 0.03 mag per the 



in-band repeatability given in Wright et al. (2010), unless the sources were brighter than 



the limits Wl = 6, W2 = 6, W3 = 4, and WA = 0; these objects were assumed to have 
magnitude errors of 0.2 mag in their respective bands. The error bar for each object's model 
magnitude was equal to the weighted standard deviation of all the Monte Carlo trial values. 
For objects with fixed t], errors on derived parameters were computed by varying rj by 0.5; 
this is approximate width of the Gaussian that was fitted to the beaming parameter for 
objects with fitted beaming (Figure [?]). For objects for which pm could not be fitted, the 
Monte Carlo trials varied Pm/pv by 1.0. 



As described in MllB and Mainzer et al. (2011c), the minimum diameter error that can 



be achieved using WISE observations is ~ 10%, and the minimum relative albedo error is 
~ 20% for objects with more than one WISE thermal band for which rj can be fitted. Table 
1 gives the results of the thermal model fits for the 428 NEOs detected by WMOPS during 
the fully cryogenic NEOWISE mission using the WISE First Pass data processing pipeline. 
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3.2. Interesting NEOs 

3.2.1. Large Amplitude Brightness Variations 

A number of NEOs stood out as objects of interest. As discussed above, we liave 
identified a number of objects with large-amplitude peak-to-peak variations. Figure [9] shows 
examples of NEOs with very large amplitude magnitude variations, indicating that these 
objects are likely to be highly elongated. WISE preferentially observed NEOs at higher 
phase angles, which may skew the observed amplitude. Effective diameters derived from 
the application of a spherical NEATM model to such objects must be viewed with caution. 
These objects would benefit both from obtaining additional observations in visible light to 
better constrain their shapes and rotational states as well as modeling them as non-spherical 
shapes. 

Many NEOs in the NEOWISE dataset are observed to have larger peak-to-peak ampli- 
tude variations in bands Wl and W2 as compared with Wi and WA. This is possibly due 
to the fact that thermal emission can emanate from beyond the terminator on an asteroid, 
depending on its thermal inertia, whereas reflected sunlight tends to come from only those 
portions that are illuminated. Since bands Wl and W2 consist of a mix of reflected sunlight 
and thermal emission, and since some NEOs have higher thermal inertias than others, a dis- 
crepancy between amplitudes at the different wavelengths can occur. An example of such an 



object can be found in NEO (1865) Cerberus (Figure 10). In addition to having significantly 
higher peak-to-peak lightcurve amplitudes in bands Wl and W2 than in bands W^ and WA, 
this object has a high beaming parameter, rj = 2.94 ± 0.03. This is close to the maximum rj 
of TT allowed by the Fast Rotating Model, which postulates that the object's cooling time is 
slow compared to its rotational rate; however, as shown in Figure [7] and as discussed above, 
rj is correlated with heliocentric distance for WISE-observed NEOs, and (1865) was located 
1.09 AU from the Sun when it was detected by NEOWISE. The object has a known rota- 



tional period of 6.810 hours (Warner, Harris & Pravec 2009), and this appears consistent 



with the rotational rate that can be observed in the figure. It may be that for this object, 
the assumption that temperature decreases as 

T{6,(f)) = Tss[max(0, cos6' cosa + sin 6* sina cos 0)]^/^ (1) 

is inappropriate, and an improved model of temperature distribution may be desirable. 
NEATM assumes that T{6, 0) goes to zero on the night side of the asteroid, even though this 
may not be correct. A number of other NEOs in our sample have high beaming parameters, 
and these objects are candidates for having high thermal inertia and/or rapid rotational 
rates. Full exploration of the true physical meaning of the relationship between rj, heliocentric 
distance, rotational rate, and properties such as thermal inertia will be the subject of future 
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work. 



3.2.2. Multi-Epoch Observations 

Due to the WISE observing cadence, there were 55 NEOs in our current sample during 
the fully cryogenic portion of the mission that were observed at two separate epochs; three 
NEOs were observed at three epochs. Of these, 20 have diameter measurements that agree 
to within 10% between all available epochs; the remaining 35 objects' diameters agree to bet- 
ter than 10% between epochs. However, three NEOs have diameter measurements that are 
disparate at the 30% or greater level. Of these, all have W2 or W?) peak-to-peak amplitudes 
larger than 0.3 mag, indicating that they are likely to have elongated shapes. Elongated 
objects observed at different viewing angles arc likely to be poorly fit by a spherical thermal 
model. Furthermore, large disparities between observing epochs can be indicative of large 
differences between temperature on the evening and morning sides; more sophisticated ther- 
mal models will need to be developed. In order to compute a single effective diameter and 
albedo for purposes of computing the debiased population statistics, the thermal model was 
fit using observations from all epochs simultaneously. These single effective diameters were 
used as an input to the debiasing model; nevertheless, these objects should be fit with non- 
spherical thermal models. The objects with relatively low amplitude lightcurve variations 
that have large discrepancies in their fitted diameters and albedos between the two epochs 
are possibly good candidates for having significant temperature differences between their 
morning and evening sides. Objects with large morning/evening temperature differences 
could be subject to increased Yarkovsky forces. 

4. Debiasing 

In order to determine the size and albedo distributions of the NEOs, it is necessary to 
remove the effects of any systematic survey biases present in the NEOWISE sample. NEO- 
WISE carried out a "blind" search for moving objects, meaning that all moving objects were 
detected in the same way regardless of whether or not they had been previously discovered 
by another observer. This means that the NEOWISE survey can be debiased independently 
of the biases of other surveys with the following caveats: There were 12 objects that were 
designated without visible follow-up and 22 objects that appeared on the NEO Confirmation 
Page but received neither designations nor follow-up; there is an undetermined number of 
tracklets that may have been misclassified by the MPC as non-NEOs. Future work will 
attempt to assess the fraction of real NEOs that are misclassified. 
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In order to model the NEOWISE survey biases, a high-fidelity simulation of the NEO- 
WISE survey was created. As described in MllA and the WISE Explanatory Supplement, 
the WMOPS pipeline required five independent detections of a particular object for it to 
have been considered detected. Sources were extracted down to 4.5 a in the Level lb images. 
We used the Julian dates and coordinates of the centers of all pointings that were used by 
WMOPS throughout the fully cryogenic portion of the survey as well as the on-sky footprint 



of each pointing (47 x 47 arcminutes) to recreate the survey history. As described in Wright 



et al. (2010), the WISE survey includes hundreds of observations of the ecliptic poles and 



an average of 10-12 observations of most moving objects in the ecliptic plane. Furthermore, 



due to the observational cadence described in MllA and Wright et al. (2010), WMOPS is 
sensitive to objects with an apparent velocity between 0.06 to 3.2 °/day. The simulated 
survey replicated these selection criteria for identifying a moving object tracklet. 

We first determined the number of NEAs with effective diameters larger than 1 km 
by collecting together as much information as possible about all known NEAs (whether 
observed by WISE or by others) with H < 22, since that is the approximate limit at which 
an asteroid with pv= 1% could still have D = 1 km. Where possible, we used diameters 
determined from radar observations, in situ spacecraft imaging, or radiometric measurements 
from WISE, IRAS, or ground-based thermal infrared observations. With NEOWISE, we 
detected 175 NEOs with D > 1 km; when we added in all the NEOs with measured diameters 
that we could find from the literature, we found an additional 75 objects, for a total of 
250 objects with previously known effective diameters that are >1 km. We also found 
from the literature an additional 72 NEAs with taxonomic classifications and assigned them 
approximate albedos based on the average albedos given in Table 1 of MllD. We note that 
most of the objects studied in MllD were Main Belt asteroids, and albedos associated with 
a particular taxonomic class were found to be strongly affected by selection biases below ~30 
km. The extent to which the albedo distributions for various taxonomic classes for the NEOs 
resemble those of large Main Belt asteroids is unknown, since the spectroscopic samples used 
in MllD and MllE are heavily biased against optically faint objects, which preferentially 
tend to be small, low albedo objects. Using the albedo distributions for different taxonomic 
classes derived from large Main Belt asteroids introduces an unknown error, and the error 
estimates we quote here must be regarded as lower limits. Our estimates will be improved 
by finding more NEOs within the WISE dataset and directly computing diameters for them. 
We used the albedo distribution of the 250 objects with known diameters to compute D for 
the previously known NEOs with H < 22 mag but no diameter or albedo measurements 
using the relationship 

'1329 ■ 10-°-2^ 
D 



Pv 



(2) 
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Applying the albedo distribution of the H <22 objects with known diameters gave us a total 
of 911±17 previously known NEOs with diameters >1 km. 

Next, it was necessary to determine the remaining population of 1 km NEOs likely to 
exist that have not yet been discovered. NEOWISE discovered 16 new NEOs with D > 1 
km. The total undiscovered population remaining was found by determining the NEOWISE 
survey bias for objects in this size range, then debiasing the NEOWISE disocoveries of large 
NEOs. The survey bias was found by creating a synthetic population of large NEOs and 
running a simulated WMOPS survey on this population to identify what would have been 
found. In order to produce a statistically valid sample, we created 25 populations each with 
50,000 synthetic NEOs with orbital elements generated randomly using the synthetic solar 



system model (S3M; Grav et al. 2011). In the S3M model, the NEO orbital elements are 



generated according to the distribution created from five source regions given in Bottke et 



al. (2002). Visible albedo and piR were randomly assigned values between zero and one; 
7] was randomly assigned a value between zero and vr. Diameter was assigned a random 
value between 1-10 km. The next step was to compute the survey bias by integrating the 
synthetic objects' orbits over the entire list of ~700,000 pointings carried out during the 
cryogenic portion of the WISE survey and see which objects were present in each frame. 
We then determined whether or not the objects would have been detected by NEOWISE by 
modeling their fluxes and comparing them to a model of the WISE sensitivity to moving 
objects in all four bands. Several other trial diameter, r], and pv distributions were used, 
and as expected, the survey bias is not sensitive to the input distributions used as long as 
there are enough synthetic objects in each bin to reduce statistical error from the simulation. 

We used the WISE Known Solar System Object Possible Association List (KSSOPAL) 
to assess the survey detection efficiency to moving objects at various places across the sky. 
KSSOPAL uses a list of known minor planet ephemerides to predict where asteroids should 
be in each WISE frame and to generate a list of probable matches; however, unlike WMOPS, 
it makes no attempt to eliminate matches to inertially fixed sources such as stars or galaxies, 
nor does it remove spurious associations with artifacts or cosmic rays. We searched the 
KSSOPAL for numbered asteroids only as these generally have well-determined orbits. In 
order to reduce the possibility of spurious associations with stars and galaxies, we checked 
each source location from KSSOPAL against the WISE Level 3 Atlas source table and used 
the n out of m statistics provided to search for sources that repeated; these sources were 
fiagged. For each magnitude bin, we computed the total number of available detections 
predicted by KSSOPAL and compared this to the total number of matches found. This 



result, shown in Figure [TT| gives an estimate of the single image completeness as a function of 
fiux for a particular region of the sky for bands W?> and WA. We computed this completeness 
curve for a number of different locations throughout the sky to bracket the WISE survey 
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sensitivity as a function of ecliptic latitude/longitude and distance from the Galactic Center. 



Figure 11 shows the probability that a moving object of a particular flux was detected by 
the WISE pipeline at a particular location on the sky. Based on the results from KSSOPAL, 
we found that the detection probability dropped to zero within five degrees of the Galactic 
Center. The detection probability curves P were fitted with the following function for both 
bands VTS and WA: 

P = ^[l-tanh(a2(W^-ai))] (3) 

where W is the W?) or WA magnitude and are the fitted coefficients. For sources between 
5 — 25° of the Galactic Center, the detection probability coefficients were given by Oq = 0.9, 
Oi = 9.5, a2 = 1.0 for band W3 and = 0.9, ai = 7.0, 02 = 2.0 for band WA; for sources 
at all further distances from the Galactic Center, the coefficients were given by = 0.9, 
ai = 10.25, 02 = 2.5 for band W3 and oq = 0.93, ai = 7.5, 02 = 2.1 for band W4. 

The synthetic population of near-Earth asteroids with D > 1 km was compared to the 
list of frames taken throughout the cryogenic portion of the survey to determine whether 
or not a particular object was present in the WISE field of view. Known comets were not 
considered in the population at this point. As discussed above, physical parameters were 
randomly assigned to each object and fluxes were computed using the NEATM model. In 
order for an object to be detected by the synthetic survey, it had to meet the WMOPS 
selection criteria: the object had to appear flve or more times and had to have a relative 
velocity in the range 0.06 to 3.2 °/day. We note that the WMOPS software is not sensitive to 
the fastest-moving NEOs, which will tend to be small objects that are very close to Earth. 
This bias against fast-moving NEOs will result in an an additional loss of sensitivity to 
the smallest objects; we will attempt to model this bias in future work. The survey bias 
for objects with D > 1 km was found by dividing the population of synthetic objects in 
each size bin by the population found by the simulated NEOWISE survey. The survey bias 



in diameter for objects >1 km is shown in Figure [13] along with the NEOWISE-discovered 
NEAs. The total number of remaining undiscovered NEAs >1 km was computed by dividing 
the NEOWISE-discovered population's differential size distribution by the survey bias, then 
integrating the result over diameter, resulting in ~80±9 NEAs remaining to be discovered 
in this size range. When added to the 911±17 previously known NEAs, we compute a total 
of 981±19 NEAs with D >1 km. This result suggests that the Spaceguard goal of detecting 



90% of all NEOs larger than 1 km (Morrison 1992) has essentially been met (although we 



note that our analysis does not yet address the near-Earth comets). 

We do not flnd strong evidence for a correlation between diameter and pv from the 



NEOWISE dataset. Harris (2006) and Delbo et al. (2003) show an apparent correlation 



between size, py, and asteroid taxonomic type. In MllD, we show that py is correlated 
with diameter within a particular taxonomic type, but this correlation could be partially or 
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entirely due to the observational biases of visible light surveys and spectroscopic measure- 
ments, which are less likely to detect and characterize small asteroids with lower albedos. 



Figure 14 shows a plot of diameters vs. py for all the NEOs observed during the cryogenic 
portion of the NEOWISE survey. No strong correlation can be observed. Furthermore, it 
should be noted that there were 12 NEOs that received designations from the MFC but 
did not have visible light follow-up, so pv could not be computed for these objects. They 
tended to be the most difficult for visible light observers to follow-up, implying that they are 
likely to have low albedos. An additional 22 candidate NEOs appeared on the NEOCF but 
were not designated due to their short observational arcs, so neither diameters nor pv could 
be computed for these objects. Like the 12 that were designated, these objects were the 
most difficult for observers to detect (however, it is possible that some of these undesignated 
objects are Main Belt asteroids or Trojans instead of NEOs). Finally, as discussed above, 
WMOPS is less sensitive to the fastest moving NEOs, which are likely to be small objects 
very close to the Earth, resulting in a further bias against the smallest objects. These fast- 
moving small NEOs are also likely to be subject to additional bias if they have low albedos, 
since we require visible light follow-up to obtain an albedo measurement. The 12 designated 
NEOs for which py could not be determined represent only ~5% of the 460 objects detected 
during the cryogenic portion of the survey, so they can only have a commensurately small 
effect on the error in the population statistics. A more detailed treatment of these objects, 
the 22 without designations, and the bias against fast-moving NEOs will be dealt with in a 
future work. For now, the error estimates presented herein must be regarded as a lower limit, 
and with these caveats, we conclude that there is no strong dependence of pv on diameter. 



The Pv distribution shown in Figure 12 was assumed to hold down to the smallest sizes. 



While a rigorous comparison of our results with those of earlier studies (Delbo et al. 2003 



Harris[2006 e.g.) would require data for NEOs of the same taxonomic types in the same size 
range, our result suggests that the apparent correlations between size and albedo reported 
by others could be entirely due to observational biases and small sample sizes as was also 



shown in Mainzer et al. (2011d). Linkages between albedo, size and space weathering cannot 



be reliably made without accounting for these biases. 

With the total number of 1 km NEAs in hand as well as the observed distributions for 
VVi Pir/pv, and r], it was possible to compute the total numbers of objects at smaller size 
ranges. As before, the total number of objects with diameters >100 m was determined by 
computing the survey diameter bias in each size bin and dividing it into the NEOWISE- 
observed diameter distribution. As was done for the >1 km NEAs, the survey biases down 
to 100 m were determined by creating a synthetic population of objects and computing 
which synthetic objects would have been detected by WMOPS. Instead of creating synthetic 
objects with flat distributions for pv, Pm/pv, and t], we used the distributions for the objects 
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that were observed by the NEOWISE survey (Figures [T2||8| and|7|). The visible albedo was 
assigned using a probability function described by a double Gaussian: 



P{p 



V 



voe 



(4) 



where Vo = 12.63, Vi = 0.034, V2 = 0.014, = 3.99, = 0.151, and = 0.122. The 
beaming parameter rj was chosen using a double Gaussian probability distribution, and 
Vir/Vv was generated using a single Gaussian probability distribution. 



D- 



The cumulative diameter distribution was initially modeled as a broken power law (A^ > 
with a = 5 above 5 km and a = 2.1 for D < 5 km based on Jedicke et al.| (2002). 



Fluxes were generated for each synthetic object using the assigned physical parameters and 
orbits as inputs to the faceted NEATM model described above. The broken power law was 
found to be a poor match to the observed cumulative diameter distribution, producing too 
many NBAs at the smallest sizes. The cumulative size distribution was found to be best 
represented by a triple power law with breaks at 5 km and 1.5 km. With this functional 
form of the cumulative size distribution, we found the slope of the power law below 1.5 
km by computing the total number of NEAs with diameters >100 m as described above 
by dividing the observed differential size distribution by the survey bias, then integrating 
the result over diameter. Using this method, we computed a total of 20,500±3000 NEAs 
with D >100 m. Using this number and the number of 1 km NEAs results in a slope of 
the cumulative size distribution power law of 1.32±0.14. This power law slope produces 
~ 13,200±1,900 NEAs with D >140 m, the limit of the George E. Brown Congressional 
mandate, and ~ 51,300 ± 11,500 NEAs with D >50 m. This result is lower than was 



predicted by previous analyses such as Harris (2008), Jedicke et al. (2002) and Rabinowitz 



et al. (2000). The break in the power law at 1.5 km is similar to that shown in Harris 



(2008). Harris (2008) shows another break in the cumulative size distribution at D ~50-100 



m and finds that the slope of the size distribution steepens, indicating that there are larger 
numbers of smaller objects in this size range. However, as our current sample only includes 
four NEAs with D <100 m, our determination of the numbers of objects or the slope of the 
size distribution below 50-100 m is not reliable. 

The survey biases in pv, Pm/pv, semi-major axis, eccentricity, inclination, and per- 
ihelion were determined by dividing the total synthetic population into the distributions of 



the objects that would have been found by the NEOWISE survey (Figure 16). The results 
show that the NEOWISE survey is essentially unbiased with respect to pv] this is expected 
given the weak dependence of thermal flux on py. Because the WISE survey covered the 
entire sky, including the ecliptic poles, the survey has a moderate bias in favor of high in- 
clination objects. Comparison of the synthetic objects' orbital elements produced by the 



Bottke et al. (2002) model to the observed orbital element distributions reveals that this 
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model is reasonably consistent with the observed population (Figure 15). Future work will 
compare the characteristics of the NEOs to their probable source regions within the Main 
Belt and comets and will attempt to further refine the distribution of NEO orbital elements. 

We finally note that it is possible that some short-arc NEOs submitted to the Minor 
Planet Center can be classified as Main Belt asteroids. It is possible to ascertain the fraction 
of objects that are misclassified by generating and sending a set of synthetic NEOs to the 
MFC in the form of tracklets (right ascensions, declinations, times, and spacecraft velocities). 
The digest score used to determine the probability that an object is an NEO could then be 
run on each synthetic tracklet to determine the fraction that would be erroneously classified 
as MBAs. This analysis will be the subject of future work in order to determine the fraction 
of misclassified NEOs likely to exist. At present, this fraction remains unquantified, and as 
stated above, the error estimates provided should be regarded as lower limits. 



4.1. NEOs As Future Targets for Exploration 



NEOs have been discussed as possible destinations for human exploration (Abell et al. 



2009), and there has been a recent push to obtain characterization information for objects 



with low Av ( Shoemaker fc Helin|1978" ), since these tend to require the least energy to reach. 



Objects with low Av are easier to reach with robotic missions as well. Figure 18 shows the 
albedo and size distributions for 80 NEOs that were observed by NEO WISE during the fully 
cryogenic portion of the mission that have Av < 7 km/s as computed by Lance Benner's list 
of Av for all NEOs observational arcs longer than 0.1 years (this being a rough measure for 
how well an object's orbit is likely to be known; objects with observational arcs shorter than 
this will most likely be extremely difficult to locate in the future). Using a similar approach. 



Mueller et al. (2011) studied 65 NEOs with Av <7 km/s; however, while there appear to be 
many targets, in practice, only a handful of NEOs will actually be suitably located in the 
timeframe of interest for human exploration (between ~2020-2050). Among our detections of 
low Av NEOs, (3361) Orpheus and (207945) have close approaches in the next two decades; 
both have moderate to high albedos (0.28 ± 0.09 and 0.16 ± 0.02, respectively) consistent 
with other S-complex objects observed by NEOWISE (MllD). NEOs with low albedos are 
potentially targets of greater interest, as these may possess volatile materials such as water 
that could be used as in situ resources for explorers. While one object among our sample 
(1996 GQ) has pv = 0.02 ±0.002 and Av ~ 6.5 km/s, it may not approach Earth sufficiently 
closely to be readily accessible to human explorers in the 2025 timeframe. 



Many of the NEOs being considered for human exploration have very large H values, 
ranging as high as if = 27 — 28 magnitudes (Abell et al. 2009). As we have shown in 



- 21 - 



earlier sections, NEO albedos range widely, from ~0.01 to ~0.6 or higher. If only visible 
photometry is available, a proposed target with H = 27 could be as small as only 5 m. 
Such an object would be similar in size (or smaller) than a visiting crew capsule: a less- 
than-ideal target. It is important to provid solid diameter estimates for any potential target 
rather than relying solely upon H due their widely varying albedos. By virtue of observing 
near 90° solar elongation, NEOWISE was able to discover an unusual NEO during the post- 
cryogenic mission that may represent the first of a new class of objects well-suited for human 
exploration. 2010 TK7 (Connors et al. 2011) is the first known Earth Trojan asteroid; 
although its inclination is too high to make its relative velocity low enough to be easily 
accessible, similar objects in more energetically favorable Trojan orbits may exist. 



5. Conclusions 

The NEOWISE project has resulted in the acquisition of a sample of NEOs that is 
essentially unbiased with respect to visible albedo, allowing us to compute the numbers, sizes, 
and albedos of NEAs down to ~100 m with reduced errors relative to previous work. By 
virtue of being in space, NEOWISE is characterized by its well-known sensitivity, number and 
location of pointings, and consistent image quality. We have shown that the Spaceguard goal 
of identifying 90% of all NEAs larger than 1 km has been met and exceeded. Furthermore, 
the number of NEAs with diameters as small as 100 m is likely to be less than previously 
suspected. We must apply similar analysis techniques to the subset of the NEAs that are 
potentially hazardous to determine whether this implies that the hazard is commensurately 
lower; this will be the subject of future work. We note that this analysis applies only to 
near-Earth asteroids; a future work will assess the population of near-Earth comets, as their 
variable levels of activity can complicate efforts to obtain sizes and albedos for them (c.f. 
Bauer et al.||2011a||bl ). 



Using the methods described above, we find that of the ~8,000 NEOs known to date, 
~5200 of them are larger than 100 m. Since we estimate that there are 20,500±3000 NEAs 
in total with D >100 m, the ~5,200 that have been found to date represent a relatively small 
fraction of the total that exist, implying that many still remain to be discovered. However, 
our sample contains only handful of objects smaller than 100 m, so we are unable to confirm 
whether or not the cumulative size distribution slope remains the same below this limit. 



Although Harris (2008) found that the size distribution slope breaks at D ~50-100 m and 
has a steeper slope below this limit, suggesting that the relative number of smaller objects 
is larger, we are unable to comment reliably on the numbers of NEAs below this size range 
with the existing NEOWISE sample at present. 
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Contrary to previous analyses with smaller datasets, we find no strong correlation of 
size with visible albedo. This result suggests that previous work (including those analyses 
attempting to link size and albedo to space weathering) has been hampered by observational 
biases, although the current lack of taxonomic classifications of most of the small NEOs 
observed by NEOWISE precludes a firm conclusion on this point (for a more extensive 



discussion of the links between size, albedo and taxonomic classifications, see also Mainzer 



et al. 2011d). We note that the inclusion of the relatively small fraction of NEOWISE- 
discovered NEOs that did not receive visible follow-up is likely to increase the numbers of 
small, low albedo objects, further reducing the likelihood of a correlation between albedo 
and size. Compared with the average distribution of visible albedos found in the Main 
Belt ( Masiero et al.||2011 ), the NEOs are preferentially brighter, following a roughly bimodal 



distribution with ~40% having pv'<0.1. As expected given that we have shown that WISE is 
essentially unbiased with respect to py, NEOWISE preferentially discovered a higher fraction 
of low albedo NEOs than that found by visible light surveys, with 53% having pv'<0.1. With 
the improved number, size and albedo distributions of the NEAs in hand, we can now begin 
to refine our understanding of their probable source regions within the Main Belt and the 
comets. 

Our thermal models have revealed that the NEATM beaming parameter 77 is correlated 
with phase angle/heliocentric distance, resulting from the fact that the temperature distri- 
bution as a function of latitude and longitude across an asteroid's surface is more complex 
than that assumed in the NEATM. Nevertheless, we have identified a number of NEOs with 
unusual properties such as (1865) Cerberus and NEOs with large-amplitude lightcurves; our 
understanding of these objects would benefit from continued study and additional follow-up 
observations. Similarly, the assumption that all objects can be reasonably represented as 
spherical should be revisited, and attempts should be made to determine three-dimensional 
shape models and rotational states. 

We note that our results will be improved by the incorporation of the second generation 
of WISE data processing, as well by continued visible light follow-up of NEOWISE-discovered 
NEOs, both to improve their orbits and to refine their H and G values. The albedos that 
we have computed depend entirely on the quahty of the H and G values that underpin 
them; while we have accounted for typical random errors associated with H and G, if these 
parameters have systematic offsets or trends, these effects can change the computed albedos. 
Our understanding of the numbers, sizes, and albedos of the NEOs will be improved by 
detection of more objects in the WISE dataset at lower signal-to-noise values as well as by 
mining the WISE data in order to compute diameters and albedos for previously known 
objects that can only be detected by combining all available images of them. We must 
also determine the debiased population statistics of the near-Earth comets and potentially 
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hazardous objects in order to better constrain the characteristics of these populations; this 
will be the subject of future work. The fraction of NEOs that are not identified correctly by 
the MPC digest score remains unquantified at present. Nevertheless, even with these caveats, 
the NEOWISE portion of the WISE project has significantly improved our understanding 
of the NEAs and paves the way for future studies. 
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Fig. 1. — WISE always surveys near 90° solar elongation; this figure shows a top-down 
view of the objects detected by NEOWISE as of February 2011 (distances are given in AU). 
The outermost circle represents Jupiter's orbit; the interior circles represent the terrestrial 
planets. Previously known NEOs are shown as green circles; new NEOs discovered by 
NEOWISE are shown as red circles; previously known comets observed by WISE are shown 
as cyan squares, and comets discovered by NEOWISE are shown as yellow squares. All other 
objects are shown as black points. The drop in density of objects observed near (+2, +2) 
AU in the figure is due to the exhaustion of the secondary tank's cryogen on 5 August, 2010, 
resulting in the loss of band WA. The dashed red line indicates the survey scan plane at the 
time of the exhaustion of the primary tank and the start of the NEOWISE Post-Cryogenic 
Mission on 1 October 2010; the survey was completed on 1 February 2011. The drop in 
detections near (+2, -2) AU in the figure is due to the intersection of the galactic plane with 
the ecliptic plane; the higher backgrounds and confusion caused by galactic cirrus resulted 
in the identification of fewer sources. 
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Fig. 2. — Follow-up observations within ~2 weeks of the last observations of NEOs candi- 
dates detected by NEOWISE were essential for securing an orbit, as the average NEOWISE 
arc spanned ~36 hours. After 2 weeks, most objects detected with such arcs have sufficiently 
high astrometric error that they cannot be recovered. This figure shows the number of follow- 
up observations contributed within 15 days of a newly discovered NEO's last observation by 
NEOWISE as a function of MFC observatory code. Spacewatch (codes 291, 691, 695 and 
807), H21 (the Killer Asteroid Froject; http://killerasteroidproject.org/student_obs.htm) , 568 
(Mauna Kea, observers D. Tholen and students) and HOI (Magdalena Ridge Observatory; 



Ryan & Ryan 2008) were the most prolific contributors within the first 15 days; other 
observers contributed follow-up over longer time spans. Observatories in the southern hemi- 
sphere are shaded black, illustrating the paucity of southern facilities that were available for 
follow-up of NEOWISE discoveries. In spite of all of the follow-up observations shown in 
this figure, 15 NEOWISE-discovered NEOs were designated without any visible follow-up, 
and another ~15-20 objects appeared on the MFC Confirmation Fage as probable NEOs 
but received neither designations nor follow-up; these objects were lost. 
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Fig. 3. — Left: WISE colors for objects observed in bands Wl, W2, W3 and WA. Black 
dots are NEOs, cyan dots are Main Belt asteroids, and green dots are Trojan asteroids. 
Right: WISE colors for NEOs observed in bands Wl, W2, W3 and WA. 




Fig. 4. — Left: W3 — WA color vs. average sky plane velocity for NEOs detected by 
NEOWISE. Right: Histogram of velocities for NEOs detected by NEOWISE. Black dots are 
NEOs, cyan dots are Main Belt asteroids, and green dots are Trojan asteroids. 
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Fig. 5. — WISE colors vs. heliocentric distance. Black dots are NEOs, cyan dots are Main 
Belt asteroids, and green dots are Trojan asteroids. 
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Fig. 6. — Most NEOs detected by NEOWISE during the cryogenic portion of the mission 
have W3 peak-to-peak amplitudes equal to ~0.4 magnitudes; however, ~18 NEOs have W3 
peak-to-peak amplitudes > 1. These objects are likely to be highly elongated and possibly 
binary. 
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Fig. 7. — Top Left: Histogram of beaming parameter [i]) values for the NEOs detected by 
the NEOWISE WMOPS pipehne during the first-pass processing carried out throughout the 
cryogenic portion of the mission. The mean value of rj for objects with fitted 1] was 1.4±0.5. 
Top Right: The beaming parameter ry vs. phase angle for cases in which the beaming 
parameter was actively fit. Black dots are NEOs, cyan dots are Main Belt asteroids, and 
green dots are Trojan asteroids. The weighted least squares fit line to the WISE NEO 
observations (red solid line) is shown compared with that of Wolters et al. (2008) (blue 
dashed line). Middle Left: The beaming parameter 77 vs. heliocentric distance. A correlation 
between rj and heliocentric distance can be observed, as shown by the red line that represents 
a running median of all the objects including NEAs, Mars crossers, and MBAs. Middle Right: 
7] VS. subsolar temperature. Bottom Left: rj vs. pv- Bottom Right: 77 vs. diameter. 
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Fig. 8. — Left: As discussed in Mainzer et al. (2011d), pm/pv can be used to distinguish 
objects with different taxonomic classifications that otherwise have similar py\ objects with 
blue or neutral visible and near-infrared (VNIR) spectral slopes such as C and B types tend 
to have lower pm/pv values than objects with red VNIR slopes, such as T and D types. 
NBAs with low pm/pv and low pv are likely to be C or B types, and objects with high 
Pir/Pv but low Pv are more likely to be D types. The shaded regions represent the regions 



covered by various spectroscopically classified Main Belt asteroids shown in Mainzer et al. 



(2011d). NEOs with known S-type spectroscopic classifications are shown as large yellow 



stars. Right: Histogram of pm/pv values for NEOs observed during the cryogenic portion 
of the mission. 
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Fig. 9. — Left: An example of an NEO with a large peak-to-peak lightcurve variation, 
(5645). This object has a peak-to-peak amplitude of 1.29 ±0.08, 1.65 ±0.06, and 1.87 ±0.08 
magnitudes in bands W2, W3 and W4, respectively. (5645) is known to have a 30.39 hour 



rotational period measured from visible photometry and is likely to be tumbling (Warner, 



Harris & Pravec 2009). A future work will convolve the visible lightcurves of this and other 



objects with the NEOWISE infrared detections. Right: NEO 2009 W06 also has a large- 
amplitude magnitude variation that is observed in all four WISE wavelengths. 
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Fig. 10. — The NEO (1865) Cerberus has a significantly larger lightcurve amplitude at 
the shorter WISE wavelengths than the longer ones: 1.31 ± 0.12, 1.28 ± 0.08, 0.54 ± 0.06, 
0.39 ± 0.07 in bands Wl, W2, W3, and WA respectively. Cerberus is known to have a 6.810 
hour rotational period ( Warner, Harris &: Pravec|2009 ), and its WISE-derived diameter is 1.6 
km. The high beaming parameter rj = 2.94 ±0.03 suggests that this could be an object with 
high thermal inertia; however, objects with low heliocentric distances tend to have higher rj 
values. The difference in amplitudes between shorter and longer WISE wavelengths could be 
due to the fact that thermal emission, which dominates bands W3 and W4, can come from 
a larger fraction of the object's total visible area than can reflected sunlight, particularly if 
it has high thermal inertia and therefore a more uniform temperature distribution. 
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Fig. 11. — The Known Solar System Object Possible Association List was used to assess 
the solar system object detection probability of the NEOWISE survey during the cryogenic 
portion of the mission across the sky. This figure shows an example of the detection prob- 
ability for a 40° x 20° region centered at 0° ecliptic latitude and 0° ecliptic longitude. The 
dots show n detected 1'^' available a fuuctiou of obscrvcd magnitude for individual numbered as- 
teroids, where Udetected is the number of times an object was actually detected out of Uavaiiabie 
possible times the object was in a WISE frame. The quantization observed at 90% and 
100% is due to the fact that on average, WISE observes most moving objects ~10 times, 
and bright objects are usually detected 9 or 10 times. The red circles are the medians of 
ndetected/navaiiabie for cach magnitude bin. The function given in Equation 2 was fitted to the 
red circles (dashed red line); this represents the detection probability. 
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Fig. 12. — Left: Preliminary observed undebiased differential diameter distribution for the 
NEOs observed during the cryogenic portion of the mission. The peak of the distribution 
occurs at ~300 m. Right: The albedo distribution for all NEOs observed by NEOWISE 
during the cryogenic portion of the mission are shown as the heavy black line; previously 
known NEOs are shown in green, and NEOs discovered by NEOWISE are shown in red. 
There are proportionately more low-albedo NEOs that were discovered by NEOWISE than 
in the previously known population; this is evidence of the bias that visible surveys suffer 
against low-albedo objects. 
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Fig. 13. — The NEOWISE survey bias for objects with D > 1 km (red hne with points) was 
computed by creating a synthetic population of NEOs with orbits generated using the orbital 
element distribution of Bottke et al. (2002). The number of NEOWISE-discovered NEOs 
>1 km is shown as the solid black line. The total number of remaining undiscovered NEOs 
was determined by dividing the NEOWISE-discovered NEO distribution by the survey bias, 
then integrating the result over diameter. 
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Fig. 14. — We do not observe a strong correlation of pv with diameter for the objects 
detected during the cryogenic portion of the NEOWISE survey. Missing from the plot are 
the 12 NEOs that received designations but received no visible hght follow-up, along with 
22 candidate NEOs that appeared on the NEO Confirmation Page but received neither 
designations nor visible follow-up. These objects are likely to be both small and dark. 
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Fig. 15. — This figure shows the properties of the synthetic population (hght gray hne 
and points) compared with the simulated objects detected by the simulated WMOPS sur- 
vey (black line and points). The medium gray hne represents the actual objects detected 
by NEOWISE during the fully cryogenic portion of the survey. All three lines represent 
normalized numbers of objects (hence the maximum value on the y axis is one). To first 
order, the orbital elements of the synthetic population found by the simulated WMOPS sur- 
vey (the "found" objects) match the properties of the objects detected by NEOWISE (the 
"observed" objects); this result indicates that the simulated objects' orbital elements are 
probably a good representation of the real population of NEAs. The good match between 
the "found" objects and the "observed" objects in py, rj, and Pm/pv is due to the fact that 
the simulated objects' physical properties were created to match those of the population 
detected by WMOPS. If we divide the synthetic population (light gray lines and points) into 
the "found" objects (black lines), the results are the survey biases for each of the various 



parameters, which are shown as the light gray lines in Figure 16 
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Fig. 16. — This figure shows the survey biases computed by dividing the synthetic population 
into the simulated objects found detected by the simulated WMOPS survey (gray lines). The 
black hnes show the properties of the objects detected by WMOPS. It can be seen that the 
WISE survey is essentially unbiased with respect to visible albedo. 
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Fig. 17. — The cumulative size distribution of the synthetic population, the synthetic 
objects found by the simulated WMOPS survey, and the objects detected by NEOWISE. 
This synthetic population of objects (gray dots) has a best-fit slope of 1.32±0.14 below 1.5 
km. The black dots represent typical error bars on the total number of objects at each size 
range. 
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Fig. 18. — There are 80 NEOs that were observed by NEOWISE during the fully cryogenic 
portion of the mission that have Av < 7 km/s and observational arcs >0.1 year (observational 
arc gives an indication of the quality of knowledge of an object's orbit), and they span a 
range of albedos similar to that observed in the general NEO population by NEOWISE. 
However, in practice, very few of these are suitable for human exploration, as they rarely 
have sufficiently close approaches in the timeframe of interest (2020-2050). 
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Table 1. NEATM results for the 428 NEOs detected by NEOWISE during the fully 
cryogenic portion of the WISE mission. An example table is given here; the full table is 
available electronically. This table contains the preliminary thermal fit results based on the 
First Pass version of the WISE data processing as described in the text. The NEOWISE 

project plans to produce an updated final catalog of physical properties based on the 
Second Pass processing of the WISE data using the updated version of the WISE Science 
Data System, with a goal of delivering this updated catalog to NASA's Planetary Data 
System. Error values presented here represent the statistical errors on the model fits, 
including Monte Carlo modeling of uncertainties for the WISE magnitudes, H, G, and 
beaming and pm when these two parameters cannot be fit. Two calibration papers 
(Mainzer et al. 2011b|[c ) discuss the absolute calibration of the WISE data for small Solar 
system bodies and should be consulted before comparing with data derived from other 
sources. The quoted precision for each parameter follows the object with the most 
significant figures for the error on that value in the table. H, G and albedo values of nan or 
"-9.99" indicate that the objects have not received visible light follow-up. Beaming value 
errors of nan or "-9.99" indicate that the thermal fit routine returned a maximum value of 
pi or a minimum of 0.3, so error cannot be properly determined. Readers are encouraged to 
check the WISE Explanatory Supplement (Cutri et al. 2011) for details and updates. 



Object 


H 


G 


D 




m) 




Pv 






V 




PIR 


N(W1) 


N{W2) 


N(W3) 


N(W4) 


K07X10C 


19.40 


0.15 


1.047 


± 


0.198 


0.028 


± 


0.014 


1.350 


± 


0.471 


0.045 


± 


0.022 








13 





K10HA8Z 


21.00 


0.15 


0.356 


± 


0.067 


0.055 


± 


0.032 


1.350 


± 


0.475 


0.089 


± 


0.051 








11 





85713 


15.70 


0.15 


3.484 


± 


0.789 


0.076 


± 


0.039 


1.350 


± 


0.450 


0.122 


± 


0.062 








9 





07822 


17.40 


0.15 


1.208 


± 


0.015 


0.133 


± 


0.022 


2.261 


± 


0.049 


0.299 


± 


0.219 


23 


24 


24 


24 


07822 


17.40 


0.15 


1.602 


± 


0.012 


0.075 


± 


0.014 


2.171 


± 


0.030 


0.139 


± 


0.043 


9 


10 


11 


11 


K02N16W 


18.00 


0.15 


0.846 


± 


0.009 


0.156 


± 


0.033 


2.118 


± 


0.066 


0.262 


± 


0.156 


30 


40 


40 


40 


K09V24O 


19.80 


0.15 


0.467 


± 


0.016 


0.098 


± 


0.020 


1.964 


± 


0.130 


0.157 


± 


0.167 





5 


5 


5 


K10MB2U 


20.60 


0.15 


0.599 


± 


0.022 


0.028 


± 


0.006 


1.684 


± 


0.111 


0.045 


± 


0.134 





7 


7 


7 


F4029 


16.40 


0.15 


2.231 


± 


0.050 


0.098 


± 


0.021 


2.044 


± 


0.085 


0.379 


± 


0.084 


5 


5 


5 


4 


F4029 


16.40 


0.15 


2.164 


± 


0.228 


0.104 


± 


0.028 


1.396 


± 


0.261 


0.166 


± 


0.045 








12 


10 


66251 


17.00 


0.15 


1.222 


± 


0.233 


0.187 


± 


0.094 


1.350 


± 


0.419 


0.300 


± 


0.150 








4 






